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Abstract 

Progress in the Effective Field Theory of two and three nucleon systems is sketched, 
concentrating mainly on the low energy version in which pions are integrated out as 
explicit degrees of freedom. Examples given are: the extraction of nucleon polarisabil- 
ities from deuteron Compton scattering at very low energies; the energy dependence 
of the nucleon polarisabilities; three body forces and the triton; and nd partial waves 
at momenta below the pion cut. 
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1 Foundations of Effective Nuclear Theory 



This presentation is a cartoon of the Effective Field Theory (EFT) of two and three nucleon 
systems as it emerged in the last three years, using a lot of words and figures, and a few 
cheats. For details, I refer to the bibliography of a recent review |1|], and to papers with 
Th.R. Hemmert @, G. Rupak §, J.-W. Chen, R.P. Springer and M.J. Savage g, P.F. Be- 
daque 0, |6|, and F. Gabbiani 0. I will mainly concentrate on the theory in which pions 
are integrated out as explicit degrees of freedom, but also comment on the extensions to 
include pions. Of the plenary talks covering related subjects, E. Epelbaum's contribution [7J 
investigates Weinberg's proposal to include pions in more detail for the few nucleon system, 
D.R. Phillips |J concentrates on electro-magnetic reactions on the deuteron, R. Timmer- 
mans |9j explains how chiral symmetry is seen in a modern potential and phase shift analysis, 
and M. Birse finally provides background on the renormalisation group point of view [|P0|j . 

For want of free neutron targets, one cannot naively extract fundamental iso-scalar and 
iso- vector properties of the nucleon separately from experiment. On the other hand, of all 
nuclei, the deuteron comes closest to an iso-scalar target and hence appears well suited to 
extract nucleon and - after removing proton effects - neutron properties. Still, the analysis 
is not straightforward because, however small the deuteron binding energy seems, binding 
effects are often not negligible at low energies where properties of the static nucleon are 
tested. Although the impulse approximation of treating the nucleons inside the deuteron 
as quasi-free is bound to become the better the higher the typical momentum scale of the 
process is, the improved resolution also necessitates a more detailed description of both the 
binding between and structure of the nucleons: Effects from meson exchanges and excited 
states are resolved at intermediate scales, and at even higher energies, the nucleons and 
mesons themselves dissolve into quarks and gluons. 

Nonetheless, model-independent predictions and extractions at low energies can succeed 
because nuclear physics provides a separation of scales. This observation is a cornerstone of 
the Effective Field Theory approach. 

Effective Field Theory methods are largely used in many branches of physics where a 
separation of scales exists. In low energy nuclear systems, the scales are, on one side, the 
low scales of the typical momentum of the process considered and the pion mass, and on 
the other side the higher scales associated with chiral symmetry and confinement. This 
separation of scales produces a low energy expansion, resulting in a description of strongly 
interacting particles which is systematic and rigorous. It is also model independent (mean- 
ing, independent of assumptions about the non-perturbative QCD dynamics): Given that 
QCD is the theory of strong interactions and that chiral symmetry is broken via the Gold- 
stone mechanism, Wilson's renormalisation group arguments show that there is only one 
local low energy field theory which originates from it: Chiral Perturbation Theory and its 
extension to the many- nucleon system discussed here. 
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1.1 The Lagrangean 



Three main ingredients enter the construction of an EFT: the Lagrangean, the power count- 
ing and a regularisation scheme. First, the relevant degrees of freedom are identified. In 
his original suggestion how to extend EFT methods to systems containing two or more nu- 



cleons, Weinberg Jl I] noticed that below the A production scale, only nucleons and pions 
need to be retained as the infrared relevant degrees of freedom of low energy QCD. Because 
at these scales the momenta of the nucleons are small compared to their rest mass, the 
theory becomes non-relativistic at leading order in the velocity expansion, with relativistic 
corrections systematically included at higher orders. The most general chirally invariant 
Lagrangean consists hence of contact interactions between non-relativistic nucleons, and 
between nucleons and pions, with the first terms reading 

£ NN = rt(iD + ^)N+ ^tr[(^Et)(^E)] + g A N^A-aN - (1) 



Coi^P^y (jV T P'iV) + [(A^PW^A^P^B - D) 2 N) +H.c. 



+ 



where N = (^) is the nucleon doublet of two-component spinors and P l is the projector 

onto the iso-scalar- vector channel, P%a = ^{ a ^)a{ T ^) b a- The iso- vector-scalar part of the 
NN Lagrangean introduces more constants Cj and interactions and has not been displayed 
for convenience. The field = y/T, = e 111 ^ 71 describes the pion with a decay constant 
f w = 130 MeV. P M = d^ + V^ is the chirally covariant derivative, and = — £^<9 M £) 

(Vfj, = |(£<9^ + C dfiC)) the axial (vector) pionic current. The interactions involving pions 
are severely restricted by chiral invariance. As such, the theory is an extension of Chiral 
Perturbation Theory and Heavy Baryon Chiral Perturbation Theory to the many nucleon 
system: The terms in the first line couple only pions amongst themselves, and to one 
nucleon, providing the familiar nuclear long range force. The terms of the second line 
couple two nucleons to each other and to pions via point-like interactions. Like in its 
cousins, the coefficients of the low energy Lagrangean encode all short distance physics - 
branes and strings, quarks and gluons, resonances like the A or p - as strengths of the 
point-like interactions between particles. As it is not possible yet to derive these constants 
by solving QCD e.g. on the lattice, the most practical way to determine them is by fitting 
to experiment. 



1.2 The Power Counting 

Because the Lagrangean (JJ) consists of infinitely many terms only restricted by symmetry, 
an EFT may at first sight suffer from lack of predictive power. Indeed, as the second part of 
an EFT formulation, predictive power is ensured by establishing a power counting scheme, 
i.e. a way to determine at which order in a momentum expansion different contributions will 
appear, and keeping only and all the terms up to a given order. The dimensionless, small 
parameter on which the expansion is based is the typical momentum Q of the process in 
units of the scale A NN at which the theory is expected to break down. Values for A NN and Q 
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have to be determined from comparison to experiments and are a priori unknown. Assuming 
that all contributions are of natural size, i.e. ordered by powers of Q, the systematic power 
counting ensures that the sum of all terms left out when calculating to a certain order in Q 
is smaller than the last order retained, allowing for an error estimate of the final result. 

For extremely small momenta Q <C m n , the pion does not enter as explicit degree 
of freedom describing long range forces. All its effects are absorbed into the coefficients 
^Cj, while formally m w — ► oo in ([!]). The only forces between nucleons are thus point- 
like two and more nucleon interactions with strengths ^Cj. This Effective Nuclear Theory 
with pions integrated out (ENT(jf), Ul2| ) was recently pushed to very high orders in the 
two-nucleon sector where accuracies of the order of 1% were obtained. It can be viewed 
as a systematisation of Effective Range Theory with the inclusion of relativistic and short 
distance effects traditionally left out in that approach. Because of non-analytic contributions 
from the pion cut, the breakdown scale of this theory must be of the order ^Ann ~ m n- For 
the ENT with explicit pions, we would suspect the breakdown scale to be of the order of 
Ma — M or m p , as A and p are not explicit degrees of freedom in ([I]). 

Even if calculations of nuclear properties were possible starting from the underlying QCD 
Lagrangean, EFT simplifies the problem considerably by factorising it into a long distance 
part which contains the infrared-relevant physics and is dealt with by EFT methods and 
a short distance part, subsumed into the coefficients of the Lagrangean. QCD therefore 
"only" has to provide these constants, avoiding full-scale calculations of e.g. bound state 
properties of two nuclear systems using quarks and gluons. EFT provides an answer of finite 
accuracy because higher order corrections are systematically calculable and suppressed in 
powers of Q. Hence, the power counting allows for an error estimate of the final result, 
with the natural size of all neglected terms known to be of higher order. Relativistic ef- 
fects, chiral dynamics and external currents are included systematically, and extensions to 
include e.g. parity violating effects are straightforward. Gauged interactions and exchange 
currents are unambiguous. Results obtained with EFT are easily dissected for the relative 
importance of the various terms. Because only S-matrix elements between on-shell states 
are observables, ambiguities nesting in "off-shell effects" are absent. On the other hand, be- 
cause only symmetry considerations enter the construction of the Lagrangean, EFTs are less 
restrictive as no assumption about the underlying QCD dynamics is incorporated. Hence 
the proverbial quib that "EFT parameterises our ignorance". 

In systems involving two or more nucleons, establishing a power counting is complicated 
because unnaturally large scales have to be accommodated: Given that the typical low 
energy scale in the problem should be the mass of the pion as the lightest particle emerging 
from QCD, fine tuning is required to produce the large scattering lengths in the S wave 
channels (1/a s ° = -8.3 MeV, l/a Sl = 36 MeV). Since there is a bound state in the Si 
channel with a binding energy B = 2.225 MeV and hence a typical binding momentum 
7 = yMB ~ 46 MeV well below the scale Ann & t which the theory should break down, 
it is also clear that at least some processes have to be treated non-perturbatively in order 
to accommodate the deuteron, i.e. an infinite number of diagrams has to be summed or 
equivalently, a Schrodinger equation needs to be solved. 
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For simplicity, let us first turn to ENT(^) in which the pion is integrated out. Here, 
a way to incorporate this fine tuning into the power counting was suggested by Kaplan, 



Savage and Wise and by van Kolck ||14j| . At very low momenta, contact interactions 



with several derivatives - like p 2 ^Ci - should become unimportant, and we are left only 
with the contact interactions proportional to $Cq. The dominating contribution to nucleons 
scattering in an S wave comes hence from two nucleon contact interactions and is summed 
geometrically in Fig. [I] in order to produce the shallow real bound state. 




i-O 



Figure 1: Re-summation of the contact interactions into the deuteron propagator. 

How to justify this? Dimensional analysis allows the size of any diagram to be estimated 
by scaling momenta by a factor of Q and non-relativistic kinetic energies by a factor of 
Q 2 /M. The remaining integral includes no dimensions and is taken to be of the order Q° 
and of natural size. This scaling implies the rule that nucleon propagators contribute one 
power of M/Q 2 and each loop a power of Q 5 /M. Assuming that 

MQ ' 2 M^A NN Q 2 ' 1 ' 

the diagrams contributing at leading order to the deuteron propagator are indeed an infinite 
number as shown in Fig. [l|, each one of order 1/(MQ). The deuteron propagator 

7' (3) 

Tnfe + A - V? - Mp - ie 



has the correct pole position and cut structure when one chooses 'Co(aO = fj 7^ • 

^Cq becomes dependent on an arbitrary scale fj, because of the regulator dependent, linear 
UV divergence in each of the bubble diagrams in Fig. [1]. Indeed, when choosing jj, ~ Q, 
the leading order contact interaction scales as in (0). As expected for a physical observable, 
the NN scattering amplitude becomes independent of fi, the renormalisation scale or cut-off 
chosen. All other coefficients ^Cj can be shown to be higher order, so that the scheme is 
self-consistent. Observables are independent of the cut-off chosen. 

The linear divergence of each bubble in Fig. [1] does not show in dimensional regulari- 
sation as a pole in 4 dimensions, but it does appear as a pole in 3 dimensions which we 
subtract following the Power Divergence Subtraction scheme |T3. Dimensional regularisa- 



tion is chosen to explicitly preserve the systematic power counting as well as all symmetries 
(esp. chiral invariance) at each order in every step of the calculation. At leading (LO), 
next-to-leading order (NLO) and often even N 3 LO in the two nucleon system, it also allows 
for simple, closed answers whose analytic structure is readily asserted. Power Divergence 
Subtraction moves hence a somewhat arbitrary amount of the short distance contributions 
from loops to counterterms and makes precise cancellations manifest which arise from fine 
tuning [14 . 



5 



1.3 Including Pions 



The power counting of the zero and one nucleon sector of the theory are fixed by Chiral 
Perturbation Theory and its extension to the one baryon sector. Therefore, the question to 
be posed is: How does the power counting of the contact terms ^Ci change above the pion 
cut, i.e. when the pion must be included as explicit degree of freedom? 

If pulling out pion effects does not affect the running of C too much, one surprising result 
arises: Because chiral symmetry implies a derivative coupling of the pion to the nucleon at 
leading order, the instantaneous one pion exchange scales as Q° and is smaller than the 
contact piece KSW C*o ~ ~ Q^ 1 . Pion exchange and higher derivative contact terms 
appear hence only as perturbations at higher orders. The LO contribution in this scheme 
is still given by the geometric series in Fig. [I]. In contradistinction to iterative potential 
model approaches, each higher order contribution is inserted only once. In this scheme, 
the only non-perturbative physics responsible for nuclear binding is extremely simple, and 
the more complicated pion contributions are at each order given by a finite number of 
diagrams. For example, the NLO contributions to NN scattering in the 1 So channel are 
the one instantaneous pion exchange and the two nucleon interaction with two derivatives, 
Fig. The constants are determined e.g. by demanding the correct scattering length and 
effective range. This approach is known as the "KSW" counting of ENT, ENT(KSW) ||13||. 




Figure 2: The next-to-leading order in ENT(KSW). 



If on the other hand, the power counting for ^Cq is dramatically modified when one 
includes pions, both the Co interactions and the one pion exchange might have to be iterated. 
In his paper, Weinberg |jnj therefore suggested to power count not the amplitude - as is 
usually done in EFT - but the potential, and then to solve a Schrodinger equation with 
a chiral potential, as pictorially represented in Fig. E. Epelbaum's talk [0] gives more 
details on the results obtained so far in this approach. 







Figure 3: The leading order in ENT (Weinberg) 



Both power countings are presently under investigation for consistency and convergence, 
and each has its advantages and shortcomings [ffl . Recently, Beane et al. fl5| demonstrated 



that Weinberg's power counting is not self-consistent in the spin doublet channel of NN 
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scattering and that pions are perturbative there, although convergence is slow. On the 
other hand, the strongly attractive 1/r 3 part of the pionic tensor force in the spin triplet 
channel necessitates a non-perturbative renormalisation of the one pion exchange. In that 
channel, the KSW counting is hence inconsistent, and Weinberg's proposal is chosen by 
Nature. It is not yet clear how this observation is to be extended to the power counting of 
counter terms in which two or more nucleons couple on external currents. After all, Nuclear 
Physics in the two body system is more than NN scattering. 

Although in general process dependent, the expansion parameter is found to be of the 
order of | in both approaches and in most applications, so that NLO calculations can be 
expected to be accurate to about 10%, and N 2 LO calculations to about 4%. In all cases, 
experimental agreement is within the estimated theoretical uncertainties, and in some cases, 
previously unknown counterterms could be determined. 

2 Applications 

2.1 Polarisabilities from Low Energy Deuteron Compton Scatter- 
ing 

The first example we turn to demonstrates not only how simple it is to compute processes 
involving both external gauge and exchange currents in a self-consistent way, but also how 
"effective" the power counting is to estimate theoretical uncertainties. The following calcu- 
lation is model-independent and performed in ENT(^) j3|. 

The iso-scalar, scalar electric and magnetic dipole nucleon polarisabilities parameterise 
the deformation of the nucleon in an external electro-magnetic dipole field: 

C pol = 27r(a E 2 + (3 B 2 )N^N , (4) 

where a := (a^W )/ 2 and (3 := ((3^+(3^)/2. How accurate does one have to measure 
elastic Compton scattering of photons on the deuteron to see them, and how severe are the 
binding effects which distinguish the deuteron from an iso-scalar target? Compare the 
graphs containing these contact interactions of a priori unknown strengths to the Thomson 
term which constitutes LO in accordance with the low energy theorem. A quick look on the 
left hand side of Fig. |] reveals that the former is suppressed against the latter by a factor 
co ,2 / ?f A^ N , where uj is the photon energy and the breakdown scale of ENT(^f) enters to get a 
dimensionless ratio. Thus, the lower the photon energy uj in Compton scattering, the less 
do a and (3 contribute. On the other hand, at uj ~ 7, the polarisabilities contribute about 
I ~ 10% to the Compton cross section, while higher energies introduce large theoretical 
errors since corrections from non-analytic pion contributions are not suppressed sufficiently 
strong any more. 

Results of an analytic calculation of the differential cross section to N 2 LO, i.e. to an 
accuracy of ~ 3%, at photon energies uj in the window of opportunity at 15 — 50 MeV are 
presented in Fig. f|[ The polarisabilities a and /3 enter as the only free parameters and 
carry a theoretical uncertainty of about 20%. As a feasibility study, we used data at uj^^ = 
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49 MeV to find a = 8.4±3.0(exp)±1.7(theor), /3 = 8.9±3.9(exp)±1.8(theor), each in units 
of 10~ 4 fm 3 . With the experimental constraint for the iso-scalar Baldin sum rule (ao + (3q = 
14), a = 7.2±2.1(exp)±1.6(theor), /? = 6.9=p2.1(exp)^1.6(theor). These values differ from 
the proton ones, a W = 12 « 10 The error bars from the experimental uncertainty in 
our extraction are large; furthermore, several conflicting measurements exist for the nucleon 
(and neutron) polarisabilities, see |§ for details. A more accurate result can be achieved (i) 
predominantly by better data, and (ii) by a higher order theoretical calculation including 
contributions from so far undetermined two-nucleon-two-photon operators. Although the 
scarcity of data at low energies is a big hindrance, the comparison of the calculation to 
experiment shows that it is quite appropriate to determine scalar polarisabilities from very 
low energy Compton scattering. The energy regime proposed is hence an interesting window 
of opportunity to determine nucleon polarisabilities in a model-independent way without 
having to deal with pions as explicit degrees of freedom. 

co = 47.72' MeV (oj Lab = 49 MeV) 
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Figure 4: Left: Thomson term and polarisabilities contribution in deuteron Compton scat- 
tering. Right: ENT(^f) result for the differential cross section, fitted to data. Dashed: «o 
and Po fitted independently; solid: «o and (3q constrained by the Baldin sum rule. The 
estimate of the accuracy of the calculation at N 2 LO (±3%) is indicated by the shaded area. 

Compton scattering has also been investigated in ENT(KSW) j|] and in ENT (Weinberg), 
see also ||. Surprisingly, the results all agree within the theoretical uncertainty of the 
calculations even at relatively high momenta. This seems to suggest that non-analytic pion 
contributions from meson exchange diagrams are small. However, the higher the photon 
energy, the less seem the polarisabilities extracted to be in agreement with the numbers 
obtained at very low photon energy, see e.g. |§. 

2.2 Dynamical Polarisabilities and Compton Scattering 

One has to keep in mind [^J the well known fact that polarisabilities themselves are energy 
dependent, probing the temporal response of the system (so-called "dynamical polarisabil- 
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ities"). One expects that the polarisabilities are enhanced around the pion production 
threshold uj ~ m n since on-shell pions can then be produced out of the virtual pion cloud 
around the nucleon. In addition, a resonance is expected at uj ~ 300 MeV because of 
contributions from the A(1232) as resonance state of the nucleon. At very low energies, 
these effects contribute corrections of order uj 2 /(m n , Ma — M^) 2 , but as uj — > m n , they are 
expected to be large. For the proton and nucleon magnetic dipole polarisability, they will 
be especially pronounced because it is anomalously small compared to and because the 
cancellation of dia-magnetic contributions from the pion tail with para-magnetic contribu- 
tions from the A and pion core which is observed at zero energy for the proton does not 
have to hold for the neutron, nor at finite energies. As a first estimate, Fig. [5] presents the 
predictions of LO Heavy Baryon Chiral Perturbation Theory (no dynamical As) and also 
estimates the effect of neglecting quadrupole and octupole polarisabilities in the extraction. 
The effect on (3 is sizeable. Thus, a series of Compton scattering experiments at energies up 
to uj ~ rn-x can give valuable information on the competing dia- and para-magnetic effects 
inside the nucleon. 




20 40 60 80 100 120 140 20 40 60 80 100 120 140 

co [MeV] co [MeV] 

Figure 5: Leading order HB^PT predictions of the dependence of the dynamical electric and 
magnetic dipole polarisabilities on the photon energy. While there is no visible dependence 
of cx.ei{oj) on the number of multipoles included (left figure), the prediction for /3 M i(uj) 
changes drastically when the extraction is truncated at the dipole polarisabilities (doted 
line in right figure), while including octupole polarisabilities makes no effect. 

2.3 Three Body Forces and the Triton 

Since all interactions permitted by the symmetries must be included into the Lagrangean, 
EFT dictates that in the three body sector, interactions like 

£ 3 bo dy = H (N^N) 3 (5) 

with unknown strength H are present. Turning again to ENT(^), we therefore have to ask 
at which order in the power counting they will start to contribute. Bedaque, Hammer and 
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van Kolck [16j found the surprising result that an unusual renormalisation makes the three 
body force of leading order in the triton channel. 

In order to understand this finding, let us first consider the LO diagrams which come from 
two nucleon interactions. The absence of Coulomb interactions in the nd system ensures 
that only properties of the strong interactions are probed. All graphs involving only ^Cq 
interactions are of the same order and form a double series which cannot be written down 
in closed form. Summing all "bubble-chain" sub-graphs into the deuteron propagator, one 
can however obtain the solution numerically from the integral equation pictorially shown in 
Fig. H within seconds on a personal computer. 




+ 




+ 




Figure 6: The double infinite series of LO "pinball" diagrams, some of which are shown in 
the first line, is equivalent to the solution of the Faddeev equation of the second line. 

Because one nucleon is exchanged in the intermediate state, each diagram is of the 
order of the nucleon propagator, i.e. Q~ 2 , while the first three buddy force (|5]) seems to 
be of order Q°. We are therefore tempted to assume that three body forces are at worst 
N 2 LO, i.e. of the order 10%. In the quartet channel, power counting suggests even N 4 LO 
because the Pauli principle forbids three body forces without derivatives. However, such 
naive counting was already fallacious in the two body sector because of the presence of a low 
lying bound state and of a linear divergence in each of the bubble graphs making up Fig. |l|. 
We therefore investigate more carefully the UV behaviour of the amplitude A(k,p) at half 
off-shell momenta p>fc,7. The homogeneous part of the integral equation simplifies then 
to 



4A(£) 
V3n 



p + q 
q p 



(6) 



where A(s) 



in the spin quartet channel, and 1 in the doublet. Qi is the Legendre 



Polynomial of the second kind, I the angular momentum of the partial wave investigated. 

The equation is easily solved by a Mellin transformation, An iS )(0,p) = p- 1+s a( l < s ) ( and 
indeed in most channels Sq is real so that only one solution exists which vanishes at infinite 
momentum and hence is cut-off independent. However, the fact that the kernel of (|6]) is not 
compact makes one suspicious whether this is always the case. Indeed, there exists one and 
only one partial wave in which two linearly independent solutions are found: the doublet 
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S wave (triton) channel, where s (l — 0, s = ~) 
any superposition with an arbitrary phase 5 is also a solution 



±1.0062 . . . i. Therefore in this channel, 



A l= n s= i(k -> 0,p) oc 



cos[1.00621np + 5] 



P 



(7) 



As each value of the phase provides a different boundary condition for the solution of the 
full integral equation Fig. |6|, the on-shell amplitude A(k,p = k) depends crucially on 5. This 
means that the on-shell amplitude seems sensitive to off-shell physics, and what is more, 
to a phase which stems from arbitrarily high momenta. A numerical study of the full half 
off-shell amplitude confirms these findings, Fig. [T[ This cannot be. 

Since physics must be independent of the cut-off chosen, this sensitivity of the on-shell 
amplitude on UV properties of the solution to (|) must be remedied by adding a counter 
term. And since the power counting in the two nucleon sector is fixed, a necessary and 
sufficient condition to render cut-off independent results is to promote the three body force 
(H) to LO, H(5) ~ Q~ 2 , and absorb all phase dependence into it, see Fig. |8|. How H(5) varies 
with the cut-off or phase is known analytically from (0), but its initial value is unknown. 
Therefore, one physical scale 5 must be determined experimentally. This one, new free 
parameter explains why potential models which provide an accurate description of NN 
scattering can vary significantly in their predictions of the triton binding energy B 3 and 
three body scattering length a^ 2 ^ in the triton channel, although all of them lie on a curve 
in the (B 3 ,a^ 2 ^) plane, known as the Phillips line, Fig. Determining H{5) by fixing 
the three body scattering length to its physical value, the triton binding energy is found in 
ENT(^) to be 8.0 MeV at LO, 8.8 MeV at NLO, and the phase shift in the triton channel 
is well in agreement with experiment [fTi 



TT] 1 — I I ■ I 1 1 1 1 1 — I I I I I ■ 1 1 1 — I I I I I ■ 1 1 1 — I I I I 1 1 1 1 1 — | 1" 
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Figure 7: Left: Variation of the cut-off by a factor 3 in a numerical study of the naive 
off-shell amplitude Fig. |6] in the triton channel changes the scattering length a(p = 0) 
dramatically. Right: Comparing the ENT prediction for the Phillips line with results from 
various potential models. Figures from Ref. ||16|| . 

It must again be stressed that the three body force of strength H(5) was added not out 
of phenomenological needs. It cures the arbitrariness in the off-shell and UV behaviour of 
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Figure 8: The cut-off independent Faddeev equation in the triton channel. 



the two body interactions which would otherwise contaminate the on-shell amplitude. Just 
as the off-shell behaviour of the two body amplitude Fig. ^ the strength H(5) is arbitrary, 
and only the sum of two and three body graphs is physically meaningful. 

Summing the deuteron bubbles, each graph in the upper line of Fig. || behaves like p~ 2 
in the UV. But the solution of the Faddeev equation goes like p~ 1+s ° with irrational (or 
even complex) Sq(1, s) |18[] and is hence more than the naive sum of graphs. 

The limit cycle thus encountered in the triton channel is a new renormalisation group 
phenomenon and also explains the Efimov and Thomas effects [TB[. In all other par- 



tial waves, three body forces enter only at higher orders [18"]. In contradistinction to the 



explanation given above, J. Gegelia's view on the triton is that a unique solution can be 



constructed using involved numerical methods, and that the three body force is not LO [19 



I thank him for intense and detailed discussions on that point during this conference, even 
though no complete agreement could be reached. 



2.4 Partial Waves in Neutron— Deuteron Scattering 

In the three body sector, the equations to be solved in ENT(^) and ENT(KSW) are compu- 
tationally trivial and can furthermore be improved systematically by higher order correction 
which involve only (partially analytic, partially numerical) integrations, in contradistinction 
to many- dimensional integral equations arising in other approaches. A comparative study 
between the theory with explicit, perturbative pions (ENT(KSW)) and the one with pions 
integrated out was performed || in the spin quartet S wave for momenta of up to 300 MeV 
in the centre-of-mass frame (E cm rs 70 MeV). As seen above, the two formulations are iden- 
tical at LO. Because three body forces enter only at high orders, this channel is completely 
determined by two body properties at the first few orders and no new, free parameters enter. 

The calculation with/without explicit pions to NLO/N 2 LO shows convergence: For ex- 
ample, the scattering length is a(LO) = (5.1 ± 1.5) fm, a(NLO, f) = (6.7 ± 0.7) fm, and 
a(N 2 LO,^) = (6.33±0.1) fm 0]. The experimental value is a( 4 S|,exp) = (6.35±0.02) fm. 
Comparing the correction of the scattering length as each order is added provides one with 
the familiar error estimate at N 2 LO: (|) 3 ~ 4%. The N 2 LO calculation is inside the error 
ascertained to the NLO calculation. The calculation of pionic corrections in ENT(KSW) 
shows that they are - although formally NLO - indeed much weaker. The difference to 
ENT(^) should appear for momenta larger than m n because of non-analytic contributions 
of the pion cut, but those seem to be very moderate, see Fig. ^|. This and the lack of data 
makes it difficult to assess whether the KSW power counting scheme to include pions as 
perturbative increases the range of validity over the pion-less theory. 
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Figure 9: Real parts of the quartet S and doublet D || wave phase shifts in nd scattering. 
Legend left: Dashed: LO; solid (dot-dashed) line: NLO with perturbative pions (pions 
integrated out); dotted: N 2 LO without pions. Realistic potential models: squares, crosses, 
triangles. Stars: pd phase shift analysis. 



Finally, the real and imaginary parts of the higher partial waves I = 1, ... ,4 in the 
spin quartet and doublet channel were found || in a parameter-free calculation in ENT(^f), 
see Fig. |9|. Within the range of validity, convergence is good, and the results agree with 
potential model calculations (as available) within the theoretical uncertainty. That makes 
one optimistic about carrying out higher order calculations of problematic spin observables 
like the nucleon-deuteron vector analysing power A y where EFT will differ from potential 
model calculations due to the inclusion of three-body forces. 



3 Outlook 



Many questions remain open: Which is the the power counting Nature chose in the pion- 
ful theory for the coupling of two and more nucleons to external currents? Does one there 
include pions perturbatively or non-perturbatively? At the moment, some people advocate a 
mixture. To investigate processes with pions in the initial or final state might be helpful |20 



How to extend the analysis of ENT(^) in the triton channel systematically to higher orders? 
Technically, how to regularise a Faddeev equation numerically with external currents coupled 
in a field theory? How do four body forces scale? Extend to nuclear matter! 
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